Spin qubits with electrically gated polyoxometalate molecules 
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Spin qubits offer one of the most promising routes to the implementation of quantum computers. 
Very recent results in semiconductor quantum dots show that electrically-controlled gating schemes 
are particularly well-suited for the realization of a universal set of quantum logical gates. Scalability 
to a larger number of qubits, however, remains an issue for such semiconductor quantum dots. In 
contrast, a chemical bottom-up approach allows one to produce identical units in which localized 
spins represent the qubits. Molecular magnetism has produced a wide range of systems with tailored 
properties, but molecules permitting electrical gating have been lacking. Here we propose to use 
the polyoxometalate [PMo^O^VO^] 9 ~, where two localized spins-1/2 can be coupled through 
the electrons of the central core. Via electrical manipulation of the molecular redox potential, the 
charge of the core can be changed. With this setup, two-qubit gates and qubit readout can be 
implemented. 



I. INTRODUCTION 

Quantum dots have often been termed artificial atoms 
and molecules.— Indeed, their tunability via electrical 
gates has permitted one to reach into a regime where 
only a single electron sits in each dot. A double dot 
then becomes the analogue of a hydrogen molecule. The 
Heisenberg exchange coupling in such a system can be 
manipulated by appropriately chosen electrical pulse se- 
quences, up to a point which has not yet been achieved 
for its real counterpart. 2 Recently, it has been demon- 
strated that this technique allows the realization of the 
fundamental one- and two-qubit quantum gates^^ in 
such a system*^ 7 - 

Naturally, the question arises whether the successful 
schemes and techniques of an electrical control of spins 
in quantum dots can be transfered back to electron spins 
in single molecules. In view of the recent experimen- 
tal progress in the field of molecular electronics, which 
demonstrates that magnetic transitions in atomic chains 8 
and single molecule o 9 ' 10 can be resolved in transport mea- 
surements, such a goal should in principle be achiev- 
able. In particular, the field of molecular magnetis m 11 ! 12 
has provided a plethora of systems of almost arbitrary 
magnetic functionality i 13 ' 14 i 15 i 16 i 17 Very recently phase- 
coherence times of up to 3 (is have been reported for such 
molecular nanomagnets after deuteration^ Yet, a suit- 
able system for an electrical control of molecular qubits 
has still been missing. It is this electrical control which 
is crucial for scalability and which distinguishes this pro- 
posal from earlier ones also based on molecular mag- 
netsM 

Here we propose an experimental setup which permits 
the electrical switching of the exchange interactions be- 
tween two electron spins. The system we choose for il- 



as an electron reservoir accommodating a variable num- 
ber of delocalized electrons hopping over the Mo centers, 
capped by two vanadyl groups containing two localized 
spins (cf. Fig. []} The spins on these two (VO) 2+ units 




lustrating our proposal is [PMoi 2 04o(VO)2] 
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a poly- 



oxometalate 17 which consists of a central mixed-valence 
core based on the [PM012O40] Keggin 22 unit, able to act 



FIG. 1: (a) Schematic drawing of the polyoxometalate 
[PMoi204o(VO)2] <J_ separated by a tunnelling barrier from 
a metallic substrate and contacted via a tunnel coupling V to 
a tip at a potential Vtt p . Indicated are the two localized spins 
Sl and Sr, of the V atoms in the center of the red square 
pyramids. Depending on the redox state of the molecule, the 
delocalized valence electrons of the Mo atoms in the center of 
the blue octahedra form a spin Sc or pair to a singlet. The 
O anions are located at the vertices of the polyhedra. (b) 
Ball-and-stick model of the polyoxometalate: O (grey), Mo 
(blue), V (red), and P (yellow). 



are weakly magnetically coupled via the delocalized elec- 
trons of the central core. We shall show how this mag- 
netic coupling can be switched in an all-electrical way and 
how this can be used for the implementation of a funda- 
mental two-qubit gate (providing entanglement), the so- 
called square-root-of-swap VSWAP. 4 Furthermore, we 
will detail how to use a variation of this procedure for 
the readout of the final state of the two qubits. 
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II. LOW-ENERGY MODEL OF THE 
POLYOXOMETALATE 

For the description of the low-energy states of the poly- 
oxometalate, two cases have to be distinguished: For an 
even number of electrons on the mixed-valence Keggin 
core, their spins pair antiferromagnetically to form a to- 
tal spin state. Then the system can be modelled by the 
two spins 1/2 on the vanadyl groups weakly coupled via 
an indirect exchange mechanism mediated by the core 
electrons. On the other hand, if the number of core elec- 
trons is odd, an unpaired spin 1/2 on the core remains 
and one obtains a set of three coupled spins 1/2. Re- 
stricting ourselves to two charge states differing by one 
electron, we can choose the electron number of the even 
state as reference and write the Hamiltonian in the form 



H = — J(nc) S L • S R - J C (S L + S R ) • s c 
+ (e - eV g ) n c + Un c (n c - l)/2 . 



(1) 



Here, the first term describes the indirect-exchange cou- 
pling between the left and right spins Sl and Sr with 
a coupling constant J(nc) which depends on the num- 
ber of electrons nc — measured with respect to the ref- 
erence number — on the central core of the molecule. In 
the case of an electron being localized on the core, its 
spin s c = (1/2) Xw d Ca T a<jl d Ccr > couples to the left 
and right spins with coupling constant Jq • Here, the op- 
erators dc,a (dc a ) destroy (create) and electron on the 
central core and r is the vector of the three Pauli ma- 
trices. The two last terms contain the orbital energy eo 
of the electron on the central core, which can be shifted 
due to a gate potential V g , and the molecule's charging 
energy U, respectively. The latter is assumed to be the 
largest energy scale of the problem. We consider a sit- 
uation where the central core of the molecule is tunnel- 
coupled to one or more metallic leads I permitting elec- 
trons to flow on and off the molecule with tunnelling rates 
IV In the STM setup depicted in Fig.QJa), two leads, the 
tip and the metallic surface, are present. We assume that 
an insulating layer — which still allows electrons to tunnel 
through — between the surface and the molecule leads to 
a molecule-surface coupling which is much smaller than 
the typical energy scale Jq of the molecular Hamilto- 
nian (TT]) . An applied bias voltage Vb — Vti p then leads to 
a shift of the molecular levels which to a very good ap- 
proximation is described by the linear relation V g = rjVtip 
with < r\ < Note that in a more sophisticated 
three-terminal setup, the parameters V g and Vb can be 
controlled independently . 24 ^ 5 



III. IDEAL QUANTUM-GATE OPERATION 



The %/SWAP operation is defined by |x, A) — > 
(Ix, A) + i|A, x)) /(I + i), where x, A =|, |. For electri- 
cally confined quantum dots, it can be physically im- 
plemented by an appropriately chosen gate- voltage pulse 



which turns on the Heisenberg exchange coupling be- 
tween the qubits for a specific time (see below) A In 
a molecular system, however, the exchange constants 
are fixed by the chemical structure and cannot be di- 
rectly changed. Thus, a more sophisticated quantum- 
gate scheme has to be employed. We propose to indi- 
rectly manipulate the coupling between the two qubits 
by changing the occupation nc of the central core in the 
molecule described above. This change in nc can be in- 
duced by changing the gate voltage V g , such that different 
charge sectors of the molecule become stable (cf. Fig. 
Chemically this amounts to a change in the redox poten- 
tial of the molecule. We will first discuss the basic idea 
of the quantum-gate scheme, thereby assuming that the 
electron number nc can be changed in a deterministic, 
externally controllable way. Later on, this assumption 
will be relaxed and the full tunnelling dynamics due to 
the change in V g will be included. 
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FIG. 2: Quantum-gate sequence for the VSWAP operation 
and exchange coupling constants during the corresponding 
gate phase. The gate is turned on for a time r ga t c . 

Initially, one adjusts the gate voltage V g in such a way 
that the stable configuration is given by an even number 
of electrons on the central core (nc = 0). We assume 
that the exchange coupling Jo = J(0) between the spins 
Sl and Sr will then be very small and can be disre- 
garded for times much smaller than H/Jq (see also be- 
low). The quantum-gate operation begins by changing 
the gate voltage such that the nc = 1 charge sector be- 
comes stable. Then, on a time-scale of the inverse tun- 
nelling rate, an electron enters the central core. After 
this has happened, the dynamics of the three-qubit sys- 
tem is governed by the Hamiltonian |T|) in the nc = 1 
subspace. The spin-dependent part of this Hamiltonian 
becomes 



H 1 = -{J 1 - J c )S r 



Sr-^S 2 . 
2 



(2) 



where S = Sl + Sr+sq is the total spin of the system and 
J\ = J(l). Before considering the time-evolution due to 
this Hamiltonian, we conclude our discussion of the gate 
cycle: after a specific time T gato , one switches back to 
the initial gate voltage and the excess electron tunnels 
off the central core again. After this has happened, the 
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two outer spins are, up to times much smaller than h/ Jo, 
decoupled again. 

The only non-trivial dynamics is induced by the Hamil- 
tonian @. The first term contains an effective exchange 
coupling between the two spins. Up to an irrelevant 
phase factor, its corresponding time-evolution yields the 
VSWAP gate if the gate time T ga t e fulfils^ 



= £ + 



(3) 



with n being an arbitrary integer. The second term in 
equation @ depends on the total spin of the three-spin- 
1/2 system. While this quantity is unknown, we can elim- 
inate its influence on the gate behavior if we restrict our- 
selves to stroboscopic gate times T ga t e for which the con- 
tribution exp [i(Jc/2)S 2 T ga te/ft)] to the time-evolution 
operator is trivial. Evaluating the effect of the operator 
S 2 in the eigenbasis (see also equation ([6]) below), this is 
the case for gate times given by 
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(4) 



fidelity T = Tr [p Ie Jt f )p ides _i] 
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with the arbitrary integer m > 0. It is important to note 
that relation ((4]) is independent of the spin direction of 
the additional electron. Equating conditions (j4| and ([3|), 
we obtain the following relation between J\ and Jq\ 
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FIG. 3: Average gate fidelity T as a function of the gate time 
Tgato and the ratio Ji/|Jc| for a ferromagnetic Jc > 0. The 
conditions ((3} and ([4]) are indicated by solid lines. The param- 
eters are: hT = 5|J C |, k B T = 0.001|J C |, and l/ g ,o n = -V g , a s = 
15|Jc| (measured from the charge-degeneracy point). 



IV. REALISTIC QUANTUM-GATE BEHAVIOR 

In the preceding discussion we have assumed that the 
tunnelling process of the electron is highly controlled, i.e., 
that we are able to instantaneously switch between the 
71q = and tic = 1 states. In reality, tunnelling is a prob- 
abilistic event occurring on a mean time scale of the order 
of the inverse tunnelling rate. In order to investigate to 
what extent this stochasticity affects our analytical find- 
ings, we compare them with numerical results obtained 
from a simulation of the incoherent quantum dynamics 
using the average gate fidelity T as figure of merit for the 
gate process (see App.lAlh In Fig.[3J we show this quan- 
tity as a function of the gate time and the ratio J\j\ Jc\ 
for the case of a ferromagnetic coupling Jc > (the anti- 
ferromagnetic case behaves very similarly). We focus on 
the strong-coupling regime fir > \Jc\i where tunnelling 
proceeds rapidly compared to the internal coherent dy- 
namics. As expected from our analytical considerations, 
the average gate fidelity assumes maxima whenever both 
conditions © and (0| are fulfilled. We found the gate 
error 1 — T to decrease with increasing tunnelling rate 
T. For T = 6|Jc|/?i, for instance, we obtain a fidelity 
T = 0.99 at the maximum indicated by the circle in 
Fig.H 



V. READOUT PROCESS 

We now show that in a molecular system described 
by the Hamiltonian ([1]), it is possible to readout the 
quantum number Sq of the two outer spins by measuring 
the sequential tunnelling current through the central dot. 
The current can be calculated using the Bloch-Redfield 
equation approach (see App. |A"|) . For a qualitative un- 
derstanding it is sufficient to consider the allowed transi- 
tions, i.e., non-zero matrix elements {a'\d c a \a), between 
the eigenstates of the two different charge sectors tlq = 
and «c = 1- The former are the eigenstates |<So, 5o, z ) 
of the total spin So = Sl + Sr. of the two outer spins. 
The Hamiltonian ([2]) of the molecule with an additional 
electron can be readily diagonalized by rewriting it up to 
an irrelevant constant as 

H^-^tf-^^Sl (6) 

Thus, in the basis of the simultaneous eigenstates 
\S,So,S z ) of S 2 , S§ and S z , the Hamiltonian in the 
«c = 1 subsector is already diagonal. For the matrix 
elements (5", S , S' z \dc CT |So, So,z) , we then obtain the se- 
lection rules S' = S ± 1/2, S' = S , S' z = S , z + a. In 
particular, we note that the quantum number stays 
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invariant under a complete sequential-tunnelling cycle. 
Since furthermore the energy differences determining the 
tunnelling rates and thus the current depend on the quan- 
tum number Sq, we find that the (quasi-)stationary value 
of the current can be strongly dependent on the initial 
preparation of So- 

This situation is depicted in Fig. Q] where the con- 
ductance peaks as a function of gate and bias voltage 
are shown schematically for two different initial values 
of Sq- From the figure we see that by starting initially 
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FIG. 4: Sketch of ground and excited state tunnel transitions 
for So ~ (dashed line) and So = 1 (solid line). 

at zero bias Vb = in the nc = state, we first reach, 
upon increasing Vb, the allowed ground-state transition 
and thus the onset of the sequential tunnelling current 
for the triplet So = 1. On the other hand, for a singlet 
preparation Sq — 0, sequential tunnelling will only be 
possible at a voltage which is higher by (J c + 2Ji)/e. 
At low enough tunnelling rate and temperature, i.e., 
HT,kBT < Jc + 2Ji, it is possible to distinguish these 
two cases and hence to measure the quantum number 
Sq. As discussed before, for an independent control of 
the two parameters V g and Vb in Fig. [H a three-terminal 
geometry is required. In the STM setup discussed here, 
one will only move along a line in Fig.|H which, however, 
is sufficient for the readout scheme, if the molecule is in 
the nc = charge sector at Vb = 0. 

Note that for an initial superposition with singlet and 
triplet contributions, the readout process just described 
represents a projective measurement. The first electron 
attempting to tunnel determines the final outcome, i.e., 
the long-time current. The average over many repeated 
measurements is described by the solution of the master 
equation (|A1[) . 

We remark that the singlet-triplet readout process re- 
quires a substantially smaller tunnel coupling than the 
one for the exchange-controlled quantum-gate operation. 
There are several ways how to achieve this: Obviously, 
one could increase the molecule-tip distance which leads 
to an exponentially strong change in the tunnelling rate. 
Still, doing so on a very short time-scale might be tech- 
nically challenging. Another possibility would be to use 
that the tunnel barrier — and, thus, in turn the cou- 



pling strength — depends on both the gate- and the bias- 
voltage. 25 Depending on the details of the contact, this 
might be enough to achieve the required change in the 
tunnel coupling. Alternatively, one could transfer the 
readout scheme to the so-called cotunnelling regimej&2! 
where transport proceeds via virtual tunnelling on and 
off the molecule. Then, the onset of the inelastic cotun- 
nelling current is solely temperature-smeared, even in the 
presence of a strong molecule-lead coupling. For measur- 
ing Sq one would then need to apply a magnetic field, 
which will lead to a splitting of the triplet state Sq = 1 
only. Thus, one will find only for this state an inelastic 
cotunnelling step in the conductance-voltage character- 
istics at the bias voltage corresponding to the Zeeman 
energy. 



VI. AB-INITIO MODELLING AND 
IMPLEMENTATION REQUIREMENTS 

While the preceding discussion applies to any sys- 
tem described by a Hamiltonian of the form fT]). 
we now return to the specific case of the molecule 
[PMoi204o(VO)2] (9_) . A high redox flexibility is a 
characteristical feature of polyoxometalates, where cyclic 
voltammetry experiments show one- or two-electron re- 
versible redox peaks, depending on the system and on the 
conditions. In particular, four compounds based on the 
XM012O40 (X = Si, Ge, P) Keggin structure (either alone 
or vanadyl-bicapped) have been recently found to show 
reversible two-electron electrochemical processes.— We 
thus extract the exchange coupling strengths for different 
charge states N of the molecule (see App. [Bjl . Identify- 
ing even and odd values of N with nc — and nc = 1 , 
respectively, we obtain the parameters in the Hamilto- 
nian ^ as given in Table U for < N < 6; for higher 
electronic populations Jo, Ji and Jc are of comparable 
size. 
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TABLE I: Orders of magnitude of the exchange coupling 
strengths in the Hamiltonian ([T]) and estimates for the size of 
the correction terms, the asymmetry parameter a = (Jcl — 
Jcr)/(Jcl + Jcr) and the gap, -Egap, between the low-energy 
states described by the Hamiltonian fT} and further excited 
states, for different electronic populations N. The charging 
energy U is of the order of 1 eV. 

It is worthwhile to consider how robust our predictions 
are against deviations from the ideal experimental setup. 
Specifically, (I) the orientation of the molecules may 
be difficult to control, (II) its electronic structure, and, 
specifically, its ideal symmetry, may be altered upon de- 
position on a surface, (III) the molecule might be affected 
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by counter ions in the surrounding, and (IV) higher- 
energy states not described by the Hamiltonian fl} might 
be energetically accessible. 

With respect to (I) we note that regardless of the orien- 
tation of the molecule, the change in charge state is firmly 
expected to affect the central M012O40 core and not the 
vanadyl groups, so we think orientation will only have 
an additional effect on (II). As an insulating thin layer 
separates the molecule from the metallic substrate, the 
alteration in the electronic structure of the polyoxomet- 
alate is expected to be minimal (details are discussed in 
the App. [B]). The main effect on the low-energy physics 
will be to introduce an asymmetry in the couplings be- 
tween the left and right vanadyl groups and the Keg- 
gin structure, Jcl and Jcr, respectively. In Table I, we 
give bounds for this asymmetry, which can be quanti- 
fied by the parameter a = (J C l - Jcr)/(Jcl + Jcr)- 
Furthermore, we have verified that an asymmetry of up 
to 10% does still allow one to achieve a fidelity of 99% 
in the regime J\ <C Jql, Jcr (f° r details see App. ICj) . 
Concerning (III), cations needed to compensate the neg- 
ative charge of the polyoxometalate polyanions, as well 
as solvent molecules, may be still present when the clus- 
ter is deposited. This feature may be a potentially im- 
portant addition on (II) which could be minimized us- 
ing an electrospray system to choose the mass/charge 
ratio corresponding to an isolated polyoxoanion. Such 
techniques have succesfully been applied to polyoxomet- 
alates and the peaks of the isolated anions have been 
detectedi 29 ' 30 Finally (IV) , the full diagonalization of the 
effective Hamiltonian also results in a set of excited states 
with energies E > E gap above those states described by 
the Hamiltonian (fTJ): we also include these data in Ta- 
ble I. 

Having these points in mind, we see that different 
charge states are usable. If we consider for instance 
N = 4 and 5, one would need molecule-lead coupling 
rates hT in the range of around 0.1 and 5meV, for the 
gate and readout process, respectively. For the gate pro- 
cedure, voltage pulses of a precision on the order of a 
few ps are required, which is not far from what can be 
achieved with present technology. In fact, if N = and 
N = 1 are found to be chemically stable in the experi- 
mental environment, their parameters will be even more 
suitable for our purposes by allowing 10 times longer 
pulse times and yielding J = as well as E gap ^> Jq,Jc- 



VII. CONCLUSIONS 

The present scheme defines a class of molecular sys- 
tems: those where two localized spins are coupled to a 
redox-active unit. The physically relevant figures of merit 
are included in Table HJ the chemical desiderata are sta- 
bility, facility of deposition on surfaces, and possibility 
of controlled oligomerization. Indeed, scalability of the 
present scheme requires covalent bonding and directed 
self-assembly of these logical building blocks. 



Molecules of this class can be found in sev- 
eral chemical families, besides bicapped polyoxometa- 
lates. Tetrathiofulvalene derivatives, polypyrrol, por- 
phyrines/phthalocyanines, fullcrene and single- wall nan- 
otubes are a few among a plethora of organic systems 
which can reversibly lose or gain one electron, can be 
substituted with radical groups, and have a rich and well- 
controlled chemistry. For example, theoretical calcula- 
tions on phthalocyanines point to the possibility of the 
chemical tailoring of these molecules to enable exchange 
couplings in the range required for the gate scheme de- 
scribed here 3 -^. 

While for an initial experimental realization with a sin- 
gle molecule, an STM contacting scheme as discussed 
in the present article enables the best control and is 
clearly favorable, a scalable method for a molecular 
monolayer could be based, e.g., on a crossbar architec- 
ture which already with current technology reaches very 
high densities^ 

In conclusion, we have proposed an experimental setup 
for single-molecule all-electric two-qubit gate and readout 
which is within reach of current technology. We have ex- 
emplified our scheme using a mixed-valence polyoxomet- 
alate, for which the model parameters have been calcu- 
lated using an ab initio approach. The general princi- 
ples behind our proposal also apply to different classes of 
molecular systems. The chemical design, synthesis and 
characterization of such systems should open new routes 
to molecular spin-qubit quantum computing. 
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APPENDIX A: BLOCH-REDFIELD-EQUATION 
APPROACH AND AVERAGE GATE FIDELITY 

For the description of the molecular dynamics in 
the presence of the molecule-lead coupling we employ 
a Bloch-Redfield-equation formalism. This approach, 
which is valid in the sequential-tunnelling regime, i.e., 
to lowest order in the tunnelling rates Ti, includes the 
full coherent quantum dynamics of the spin state a 33 ! 34 
as described by the density matrix elements Q a p = 
Tri eac is(Q;|p|/3) after tracing out the leads. In the usual 
Born-Markov approximation, we obtain 

Q a p = - iu a f3Q a p + - ^ \[ W f3>0aa> + { W L'af3p' )*] Qa' fi' 
la'? 

- W l/3'0'a' Qa'f) ~ {Wp a , a ,p, )* Q a p' j , 

(Al) 
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where we have introduced the frequencies LJ a p = (E a — 
Ef})/h and the transition rates Wp l/3aa , = WJg aa , + 

Wp7p aa , due to tunnelling across a molecule-lead con- 
tact £. The rates for tunnelling of an electron on and 
off the molecule assume the form Witg aa , = r> f(E a — 

E a '-M)EAP'\dcMPMdc,>') and^ QQ , = T t [1- 
f(E a , -E a - fit)} E £ r(^'MLl/ 3 X Q! Mc,a|a'), respectively. 
Here, we assume that the lead electrons are described 
by a Fermi distribution /(e — m) = {1 + exp[(e — 
/^)/£;T]} _1 at temperature T and electro-chemical po- 
tential [/,£. The tunnelling current across contact I can 
be determined from the density matrix elements via 
h = eReE Qa / (W%- Ia , a - W l + a , a , a )g aP . 

Using the Bloch-Redfield equation approach, we are 
able to simulate a realistic gate cycle. This al- 
lows us to calcu late the average gate fidelity 35 T = 
Tr [preai (*f )pideal] as quantitative measure for the gate 
quality. Here, the overline indicates the average over 16 
unentangled input states I^^lI^^r, where l^i) = ||), 
|* 2 > - 14), |*3> = (IT) + \i))/V2, and | tt 4 ) = (IT) + 
i|4))/V2, Pideai is the ideal result of the %/SWAP oper- 
ation as defined above and p re ed{tf) is the state of the 
system according to the quantum dynamics (|A1|) at the 
final time t j . 

APPENDIX B: EVALUATION OF EXCHANGE 
PARAMETERS 

For the evaluation of the exchange coupling constants 
and to check the validity of our numerical assumptions, 
we used an effective Hamiltonian approach. We consid- 
ered the main one- and two-center phenomena in mixed- 
valence systems, namely magnetic exchange, electron 
transfer, Coulomb repulsion and orbital energy, which 
were parametrized for [PMoi204o(VO)2] ? ~ (details will 
be published elsewhere). Thus, we diagonalized a 14- 
site t—J—V-e effective Hamiltonian for different electronic 
populations. We projected this full energy level scheme 
into a subsystem containing only the localized spins in 
the [V /v O] groups, plus, for odd values of N, an un- 
paired spin in the M012O40 Keggin structure. 

For N even, Jq is given by the energy difference be- 
tween the lowest singlet and triplet states, while E gap is 
the gap to the next energy level. For N odd, the evalua- 
tion of E gap is equally direct. Obtaining the exchange pa- 
rameters is also straightforward when a symmetric molec- 
ular structure is assumed: Jq and J\ can then be readily 
calculated from the energy differences between the first 
doublet and quadruplet states. When deviations from 
symmetry are considered, the determination of the three 
exchange parameters Jq , Jcl and Jcr requires the anal- 
ysis of the wavef unctions. 

The spatially-dependent microscopic parameters t, J, 
V, and e were obtained from ab-initio calculations 
of molecular subclusters, which accounts for inhomo- 
geneities in the molecule. The experimental setup is ex- 



pected to imply small variations in bond distances and 
angles, as well as a different charge distribution in the 
surrounding of the molecule. In order to assess how 
these perturbations affect the parameters in the effec- 
tive Hamiltonian the calculations were repeated for 
a reasonable range of deviations (±50% in J and t, and 
±leV in e). Except for the case with N — 3, where the 
excited states not contained in the Hamiltonian |TJ lie 
close in energy, the parameters of our interest proved to 
be robust enough within the orders of magnitude indi- 
cated in Table HI 

APPENDIX C: CHARACTERISTIC GATE 
QUANTIFIERS 

Besides the average gate fidelity discussed in the main 
text, one can consider further quantities for assessing the 
quality of a quantum-gate operation 3 ^. One is the so- 
called gate purity 

P = Tr[p roa i(t f ) 2 ], (CI) 

which is a measure for the decoherence of the gate. 
Here, the average is over all 16 unentangled input states 
IV'uo(io)) = |*<)l|*j)r, i,j = 1,...4 (seeApp.EJfor the 
definition of the l^j)). For an ideal \/SWAP operation 
it reaches unity. Another measure is the entanglement 
capability C which is defined as the smallest eigenvalue 
of the partial transpose of the density matrix^ p lea .i(t{) 
obtained for all 16 unentangled input states just defined. 
Ideally, it assumes a value of —1/2. Finally, the quantum 
degree of the gate is given by 

Q= max (V'molPreal(if)IV'mc) , (C2) 

|Vmo),Puc(t ) 

i.e., as the maximum of the overlap between all pos- 
sible output states obtained for all 16 unentangled in- 
put states and all maximally entangled states \ip m e), 
viz, the four Bell states (|T)l|T)r ± |T)l|T)r)/\/2 and 
(U)l|T)r ± |T)l|T)r)/%/2. In the optimal case, the quan- 
tum degree reaches 1/2 for the vSWAP gate. 

In Fig.O we show these four gate quantifiers as a func- 
tion of the gating time T ga t e and the ratio J\/\Jc\- The 
gate fidelity T has already been discussed in the main 
text. We observe that the gate purity V assumes almost 
unity for the stroboscopic gating times (4). This means 
that at these times, the additional spin sc and the spin 
S = Sl + Sr. become disentangled again, such that one 
recovers a pure state when removing the additional elec- 
tron. Furthermore, we find that all gate quantifiers ex- 
cept for the gate fidelity overestimate the quality of the 
gate: They almost assume their perfect values even for 
parameters where the gate fidelity is much below unity. 
The quantum degree Q turns out to be completely in- 
dependent of the gating time and the exchange-coupling 
ratio. 

The situation in the presence of an asymmetry a — 
( Jcl - Jcr)/ {Jcl + Jcr) between the couplings J C l and 
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FIG. 5: Four average gate quantifiers (a) gate fidelity (b) purity (c) entanglement capability, and (d) quantum degree each 
normalized to the maximal value achievable for an ideal VSWAP operation as a function of the gating time r ga to and the ratio 
Ji/\Jc\- The parameters are as in Fig. 3. The conditions (3) and (4) are indicated by solid lines. Note that the color scale for 
the gate fidelity is different from the one used in Fig. 3. 
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Jcr of the central Keggin core to the left and right, re- 
spectively, vanadyl groups is shown in Fig. [BJ We have 
chosen the upper bound a = 0.1 from the results of our 
ab initio calculations (cf. Table I). We find for the gating 
times shown in Fig.[6]that this noticeable aberration from 
perfect symmetry does not strongly affect the gate qual- 



ity for small and/or negative values of the ratio J\/\Jc\, 
where Jq = (Jcl + <^cr)/2 is the average coupling of the 
central spin to the left and right spins. In particular, we 
still find a fidelity better than T = 0.99 at the maxi- 
mum indicated by the circles in Fig. [6] for a tunnel rate 

r = ii|J c |/R. 
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FIG. 6: Four average gate quantifiers (a) gate fidelity (b) purity (c) entanglement capability, and (d) quantum degree each 
normalized to the maximal value achievable for an ideal VSWAP operation as a function of the gating time T gato and the ratio 
■h/\Jc\, where Jc = (Jcl + Jca)/2 is the average exchange coupling between the Keggin core and the left and right vanadyl 
groups. An asymmetry a = (Jcl — Jcr)/(Jcl + Jcr) = 0.1 of these exchange couplings is assumed. The other parameters are 
as in Fig. 3. The conditions (3) and (4) are indicated by solid lines. 



